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ABSTRACT 

The problem of scheduling sensor transmissions for the detection of corre- 
lated random fields using spatially deployed sensors is considered. Using 
the large deviations principle, a closed-form expression for the error ex- 
ponent of the miss probability is given as a function of the sensor spacing 
and signal-to-noise ratio (SNR). It is shown that the error exponent has a 
distinct characteristic: at high SNR, the error exponent is monotonically 
increasing with respect to sensor spacing, while at low SNR there is an 
optimal spacing for scheduled sensors. 

1. INTRODUCTION 

One of the most critical design constraints for large scale sensor 
networks is energy efficiency. In the context of detecting spatially 
correlated signals, this means that, given a desired level of detector 
performance, one should minimize the amount of data required, 
and collect data from judiciously chosen areas. 

We consider in this paper the detection of a one dimensional 
diffusion process that has a spatial state space structure as illus- 
trated in Fig.Q Specifically, we assume that the underlying signal 
s{x) is the stationary solution' of the diffusion equation 



ds{x) 
dx 



-As{x) + Bu{x), x>Q, 



(1) 



where A> Q and B are known, and the initial condition is given 
by s(0) ~ A/'(0, Ho). The input process u{x) is zero-mean white 
Gaussian, independent of both sensor noise and s(0). 




Fig. I. Signal and sensor location 
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' The stationary solution of Q requires a condition on A, B, Ho given 

byno = 1^. 



We assume that each sensor in the field takes a noisy measure- 
ment. At location Xi, the sensor measurement yi comes from the 
following binary hypotheses 



Ho 
Hi 



2/i = Si + Wi 



1,2, 



(2) 



where Si — s{xi), and Wi are i.i.d. sensor measurement noises 
from A/'(0, cr^) with a known variance . We will not focus here 
on minimizing the number of bits for the quantization of t/i, nor 
will we tackle the problem of how yi are sent to the fusion cen- 
ter for detection. These are important design issues that must be 
treated separately. Our focus in this paper is that, if data are to be 
collected from the sensor field, where they should be collected. For 
energy efficiency, we aim to find a scheduling scheme that require 
as few samples as possible. 

I.I. Summary of Results 

We adopt the Neyman-Pearson formulation by fixing the detector 
size a and minimizing the miss probability when the number n 
of samples is large. Specifically, we are interested in two closely 
connected design problems: the locations X = {xi} where data 
should be collected, and the rate of decay of miss probability for 
detectors of size a. We will make the assumption that samples 
are collected uniformly with equal spacing = A. Such an 
idealized assumption can only be approximated in practice, but 
it does provide an analytically tractable formulation that leads to 
insights into energy efficient data collection. 

The miss probability Pa/ (A, n; a, SNR) is a function of sen- 
sor spacing A, the sample size n as well as detector size a and 

SNRT = %. The energy efficient scheduling problem can then 
be formulated as one of optimizing the error exponent Ka (A, SNR) 
under the large deviation principle where 



K^{A,T) = lim -IogPAf(A,n;Q,r). 

n — >oo Tl 



(3) 



The connection of the above formulation with energy efficient sen- 
sor scheduling is natural when the sample size n is directly related 
with the number of transmissions (such as the case in Sensor Net- 
works with Mobile Access (SENMA)). 

We derive a closed-form expression for the error exponent 
ifct(A, r) of miss probability (which is independent of a) by ex- 
ploiting the state-space structure of alternative hypotheses and mak- 
ing a connection with Kalman filtering. We show next that Ka (A, V) 
has a distinct phase transition: when SNR F > 1, Ka{/^,T) is 
a monotonically increasing function of A, indicating that sensor 
spacing A should be made as large as possible (for fixed but large 
sample size). When SNR F < 1, on the other hand, Ka{A.,V) 
achieves the maximum at some A*, which means that there is an 
optimal spacing among collected samples. We also provide an 



implicit equation for A*. We also present simulation results to 
demonstrate the predicted behavior. 

1.2. Related Work 

The detection of Gauss-Markov process in Gaussian noise is a 
classical problem. See 1 5 1 and references therein. Our work re- 
lies on the connection between likelihood ratio and the innovation 
process through Kalman filtering by Schweppe |^. While the con- 
nection between Kalman filter and error exponent is a contribution 
of this paper, there is an extensive literature on the large devia- 
tion approach to the detection Gauss-Markov process 1 9**101 flTI 
\l2 13 14 15 1. These results do not provide explicit expressions 
(with an exception for noiseless AR processes) from which opti- 
mal scheduling can be obtained. 

The sensor scheduling problem can also be viewed as a sam- 
pling problem. To this end, Bahr and Bucklew ['11 1 optimized the 
exponent numerically under a Bayesian formulation. For a specific 
signal model (low pass signal in colored noise), they showed that 
the optimal sampling depends on SNR, which we also show in this 
paper in a different setting. 

2. SYSTEM MODEL 

We derive the discrete model for the sampled data of the signal 
s{x) similarly to [4|. For equally spaced sensors with spacing A, 
the dynamics of sampled signal are described by 



Si+i = asi + Ui 



(4) 



a 



Ui — e~ ^ Bu{a 
Jo 



(Xi+i — x)dx, 
The mean and variance of process Ui are given by 

Em, = 0, Q = Var(uO = no(l-a2). 



(5) 



The signal samples {si} form an autoregressive sequence under 
Hi. Notice that Esf — flo for all i and the value of a determines 
the amount of correlation between signal samples. For i.i.d. signal 
samples we have a = and a = 1 for the perfectly correlated 
signal. Notice also that the noisy observation {t/;} are not autore- 
gressive; they follow the hidden Markov model. 

3. ERROR EXPONENT 

In this section, we investigate the performance of Neyman-Pearson 
detector with a level a G (0, 1). We present a closed-form expres- 
sion of the error exponent of miss probability Ka{^, F) defined 
in J3}- The theorem below comes from the fact that the limit of the 
normalized log-likelihood ratio under Hq is the best error exponent 
for general ergodic cases EH. By expressing the log-likelihood 
ratio through the innovation representation, we make the closed- 
form calculation of error exponent tractable. 

Theorem 1 (Error exponent) For the Neyman-Pearson detector 
of the hypotheses J2| with level a £ (0,1) (i.e. Pf < a) 
and < a < 1, the best error exponent of miss probability is 
independent of a and is given by 



i^<,(A,r) = -iiog — + 

^ ^ 2 ^ Re 2 Re 2 



(6) 



where Re and Re are the steady-state variances of the innovation 
process ofyi calculated under H\ and Ho, respectively. Specifi- 



cally, Re and Re are given by 



P2 + 2cr2p + (1 - a2)CT4 



y [<t2(1 - a2) - Q]2 + 4a2Q - cr^(l - 



(7) 
(8) 

(9) 



In frequency domain. 



7^,(A,r) = -L r D{N{Q,cT^)\\N{Q,Sy{uj)))duj, (10) 
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where D{-\\-) is the Kullback-Leibler distance, and the spectrum 
Sy{uj) of{yi} under H\ is given by 



Sy{uj) 



2 I 



1 — 2a cos uo + 



(11) 



Proof: see 1171 . ■ 

For notational convenience, we use K for the error exponent. 

3.1. Properties of Error Exponent 

First, it is easily seen from TheoremQthat K is a continuous func- 
tion of the correlation coefficient a (0 < a < 1) for a given IIo 
and (7 since i?e > cr^ > 0. 

Theorem 2 The error exponent is positive for any SNR V and < 
a < 1. Furthermore, 

(i) for i.i.d. observations (a = 0), the error exponent reduces 

to the Kullack-Leibler distance i3(po| where po ~ A/'(0, cr^) 
andpi ~ A/'(0,no + o-^); 

( ii) for Perfectly correlated signal (a — 1), the error exponent 
is zero for any SNR F, and the miss probability is bounded 
by 



1 



27V 



D)c 



-1/2 



< Pai < 



1 



-1/2 



(12) 



for sufficiently large n, where c and D G (0, "^^) cire 
positive constants. 

The positivity of K is immediate from Theorem Q The case 
when a = Q corresponds to the Stein's lemma for the i.i.d. case. 
Unless a = 1, the miss detection probability always decays expo- 
nentially. 

For the perfectly correlated case (a = 1), the miss probability 
does not decay exponentially; it decays with 0(--^) as shown in 
<12t . This is explained by the form of the optimal detector. It can 
be shown that the sufficient statistic is given by 



(13) 



Under Hq, we have X]"=i Vi ~ •'V(0, na"^). Since f(x) = x^ is 
symmetric about zero, the Neyman-Pearson detector with level a 
is simply given by 



(14) 



where Zn = \fnaQ (Q(-) is the tail probability of stan- 

dard normal distribution). Since Pi,„ — Af{0, n^Ho + na"^), the 
probability of missing is given by 

^« = l-20 4=== • (15) 
y V w^IIo + ncr^ / 

For large n, Pm behaves as 1 — 2(5 (-^) where c is a constant, 
which leads to <12> . 

Having obtained the behavior of error exponent at two extreme 
correlation cases, we now show that the error exponent has distinct 
characteristics at different SNR regimes. 

Theorem 3 {K vs. correlation - high SNR ) K is monotone de- 
creasing as the correlation strength increases (i.e. a ] 1) if the 
SNR r > 1. 

The above theorem shows that the i.i.d. observations give the 
best error performance for the same SNR with the maximum error 
exponent D(J\f{o, l)||A/'(0, l + r)) when SNR is larger than unity. 
Fig. |2|(left) shows the error exponent as a function of the corre- 
lation coefficient a. In the context sensor scheduling, if n trans- 
missions are required for the desired miss probability, it is best to 
maximize the distances between the scheduled sensors. The intu- 
ition is that, at high SNR, the signal component in the observation 
is strong, and the innovation process contains more information 
about the hypotheses. Thus there is benefit to make the signal less 
correlated. 




Fig. 2. K vs. a (left:SNR=10 dB, right:SNR=-3,-6,-9 dB) 



In contrast, error exponent at low SNR is not monotonic, and 
there exists an optimal correlation. Fig. |2| (right) shows the error 
exponent for the case that the SNR is less than unity. As shown, 
the i.i.d. case is not the best for the error performance for the same 
SNR. The error exponent initially increases as a increases, and 
then decreases to zero as a approaches one. For SNR of -6 dB the 
error exponent is less than the case of -3 dB and it is seen that am 
is shifted closer to one. The following theorem characterizes the 
optimal correlation. 

Theorem 4 (K vs. correlation - low SNR) There exists a non-zero 
value am of the correlation coefficient that achieves the maximum 
K for SNR < 1, and am is given by solving the following equa- 
tion. 

[l + a2+r(l-a^)]^-2(re + — ) =0, (16) 

where = R^/o^. Futhermore, am converges to one as SNR V 
goes to zero. 



At low SNR, noise in the observation dominates. Thus making 
signal more correlated provides the benefit of noise averaging. The 
lower the SNR, the stronger the correlation is desired as shown in 
Fig.|2|(right). Note however that excessive correlation in the signal 
doesn't provide new information by the observation. Notice also 
that the maximum value of error exponent is much larger than the 
value for i.i.d. case at lower SNR. Hence, the optimal correlation 
gives much better performance with the same number of sensor 
observations for low SNR cases. 

A significance of the above theorem is the determination of 
optimal sensor spacing. In particular, the optimal distance is given 
by 

A* = -log(a™)/A, (17) 
for the same underlying physical phenomenon described by Q. 




Fig. 3. Left: a,„ vs F, rightiA* vs r(< 1) 

Fig. |3|shows the value of am that maximizes the error expo- 
nent as a function of SNR. As shown in the figure the SNR of unity 
gives sharp transition between two different behaviors of error ex- 
ponent w.r.t. correlation strength. 

Finally, we investigate the behavior of the error exponent w.r.t. 
SNR. 

Theorem 5 (K vs SNR) The error exponent K is monotone in- 
creasing as SNR increases for a given correlation coefficient < 
a < 1. Moreover, at high SNR K increases linearly with log SNR. 




Signal-to-Noise Ratio [dB| 



Fig. 4. K versus SNR (a = e"^) 

The log SNR increase of K w.r.t. SNR is similar to the case of 
diversity combining of Rayleigh-faded multipaths in additive noise 
since in both cases the signal component is random. Comparing 
with the detection of a deterministic signal in noise where the error 
exponent is proportional to SNR, the increase of error exponent 
w.r.t. SNR is much slower for the case of stochastic signal in noise. 

3.2. Sensor Placement 

So far we have assumed that sensors have been placed and the 
problem is to schedule the transmission of selected sensors. A 



related problem is sensor placement. Should n sensors be placed 
to cover as large an area as possible, or should they be clustered in 
a subregion of the signal field? Theorem |3| suggests that, at high 
SNR, these sensors should be placed with maximum separation 
subject to the size of the field. At low SNR, however, sensors 
should be placed with the optimal separation. This implies that, if 
the sensor field is large, it is better to cluster the sensor in some 
region, leaving other areas without sensors. 

4. SIMULATION RESULTS 

We considered the Neyman-Pearson detector of first-order autore- 
gressive signal described by ^4)- We considered SNR of 10 dB and 
- 3 dB, and several a for each SNR. The probability of false alarm 
was set 0. 1% for all cases. Fig.|5|shows the miss probability w.r.t. 




Fig. 5. Pm vs. # of sensors (SNR=10dB) 

number of sensors for 10 dB SNR. It is shown that the i.i.d. case 
(a = 0) has the largest slope for error performance, and the slope 
of error decay is monotonically decreasing as a increases to one. 
Notice that the error performance for the same number of sensors 
is significantly different for different correlation strength even for 
the same SNR, and the performance for weak correlation is not 
much different from i.i.d. case predicted by Fig. |2|(left). (We can 
see that the slope decreases suddenly near a = 1.) It is shown that 
the miss probability for highly correlated case (a — 0.9999) is not 
exponenitally decay. 




Fig. 6. Pm vs. # of sensors (SNR=-3dB) 

The error performance for SNR of -3 dB is shown in Fig. |S| 
It is seen that the slope increases as a increases from zero, and 
reaches maximum with sudden decrease after the maximum. No- 
tice that the error curve is still not a straight line for low SNR case 
due to the o(n) term in the exponent. Since the error exponent 
increases only with log SNR, the required number of sensors for 
-3 dB SNR is much larger than for 10 dB SNR for the same miss 
probability. It is clearly seen that Pm is still larger than 10^^ for 



200 sensors whereas it is 10 * with 20 sensors for 10 dB SNR 
case. 

5. CONCLUSIONS 

We have considered the detection of correlated signal using noisy 
sensors. We have derived the best error exponent for the Neyman- 
Pearson detector satisfying a fixed size constraint using the inno- 
vations and the spectral domain approaches. We have also inves- 
tigated the properties of the error exponent. The error exponent is 
a function of SNR and correlation strength. The behavior of er- 
ror exponent w.r.t. correlation strength is sharply divided to two 
regions depending on SNR. For SNR larger than unity the error 
exponent is monotone decreasing as correlation becomes strong. 
On the other hand it has a non-i.i.d. correlation strength that gives 
the maximum slope for SNR smaller than one. Using the property 
of error exponent, the optimal strategy for sensor scheduling has 
been derived. 

The results presented in this paper have a generalization to the 
vector state-space model. Such a generalization is useful to con- 
sider more general scheduling schemes. In practice, the locations 
of sensors are random, and sensors may cease to function. Thus 
random sampling may need to be considered. 
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